function kap = kap_pscc(shat,kap0,phi)
%selling cost per existing customer
kap = kap0/phi*shat.^phi;
end

